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A flux liquid can condense into a smectic crystal in a pure layered superconductors with the 
magnetic field oriented nearly parallel to the layers. If the smectic order is commensurate with 
the layering, this crystal is stable to point disorder. By tilting and adjusting the magnitude of the 
applied field, both incommensurate and tilted smectic and crystalline phases are found. We discuss 
transport near the second order smectic freezing transition, and show that permeation modes lead 
to a small non-zero resistivity and large but finite tilt modulus in the smectic crystal. 
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In the past few years, both experimental [gjj^ and the- 
oretical [H^] work has emphasized the importance of pin- 
ning in type II superconductors. Although much of this 
work has focused on random defects, e.g. twin bound- 
aries, the layered structure of the copper-oxide materi- 
als itself provides a non-random source of pinning 
At low temperatures, since the c-axis coherence length 



cO 



4A < 



12A, the lattice constant in this direc- 



tion, vortex lines oriented in the ab plane are attracted 
to the regions of low electron density between the Cu02 
layers. In this paper, we discuss the phase diagram of in- 
trinsically pinned vortices in near perfect alignment with 
the layers. In contrast to previous work, we focus on the 
behavior relatively close to Tc, where transport measure- 
ments are most easily performed and hysteretic effects 
are weak. Our research is motivated by the recent ex- 
perimental work of Kwok et. al. who observed a 
continuous resistive transition in YBa2Cu307 for fields 
very closely aligned {6 < 1°) to the ab plane. 

To understand thermal fluctuations, we employ the 
boson mapping 0. Consider an isolated flux line ori- 
ented along the a-b plane. For T > 80K, ~ '?co(l ^ 



-1/2 > 



s, and the intrinsic pinning barrier energy 



per unit length is 



[/p « 5 X 10^^ ( ^ 
7 V s 
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-15.8Cc/s 



(1) 



where the energy scale eg = {(po/ATrXab)'^ , with (j)Q the 
flux quantum and Xab the a~b plane penetration depth, 
and where 7 = ^ nic/mab is the anisotropy. With coor- 
dinates y II B _L c, z II c, the free energy of the vortex 
(described by x{y) and z{y)) is 
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where the stiffness constants obtained from anisotropic 
Ginzburg-Landau (GL) theory are e|| = eo/7 and e± — 
eo7 (see, e.g. Ref. 0). Upon integrating out the x- 
displacement, the remainder of Eqj2| maps to the Eu- 
clidean action of a quantum particle in a one-dimensional 



periodic potential. In the quantum-mechanical analogy, 
the particle tunnels between adjacent minima of the pin- 
ning potential, leading to the well-known completely de- 
localized Bloch wavefunctions even for arbitrarily strong 
pinning. The "time" required for this tunneling maps 
to the distance, Lkink, in the y-direction between kinks 
in which the vortex jumps across one Cu02 layer. The 
WKB approximation gives 



L 



kink 




Ups/k^T 



(3) 



where fceT plays the role of h. When the sample is larger 
than -Lkink along the fleld axis, the flux line will wander 
as a function of y, with 



([^(y)-^(o)]^ 



Dy, 



(4) 



where the 



'diffusion constant 

< 



D sa s^/Lkink- 
the pinning is extremely weak, 



For ^i±UpS < ksT, 
and the WKB approximation is no longer valid. Instead, 
the diffusion constant D w fceT/ej^, as obtained from 
Eq.^ with Up = 0. At much lower temperatures, when 
<C s, the energy in Eq.|l] must be replaced by the cost 
of creating a "pancake" vortex ||l^ between the Cu02 
planes. In this regime, Lkink ~ Lbis/^cY""^'''''^ ■ 

For T w 90K, as in the experiments of Kwok et. 
al. §, Cc/s ~ 2.3, and Eq|| gives yO^s/fceT < 1, 
indicative of weak pinning and highly entangled vor- 
tices in the liquid state. The transverse wandering in 
this anisotropic liquid is described by a boson "wave- 
function" with support over an elliptical region of area 
k-QTLy/ ^e\\e± with aspect ratio Ax/Az = 7^5 for 
YBa2Cu307 {Ly « 1mm is the sample dimension along 

y) 0. ^ 

To explain the observed transition, the interplay be- 
tween inter-vortex interactions Q and thermal fluctu- 
ations must be taken into account in an essential way. 
The experiments of Ref. Q rule out conventional freez- 
ing, which is flrst order in all known three-dimensional 
cases. A vortex/Bose glass (VG/BG) transition also ap- 
pears untenable, due to the purity of the sample (only six 



1 



twin planes are present, and first order melting is seen for 
> 1°). The dynamical scaling exponents are also in- 
consistent with VG/BG values. Instead, we postulate 
freezing into an intermediate "smectic" phase between 
the flux liquid (high temperature) and crystal/glass (low 
temperature). Such smectic freezing, as first discussed 
by de Gennes for the nematic-smectic A transition [ p2[ , 
can occur via a continuous transition in three dimen- 
sions. Changes in the applied magnetic field tune the 
commensurability of the vortex system with the under- 
lying periodic structure . 

Our analysis leads to the phase diagram shown in 
Fig.|l|. Upon lowering the temperature for He = and a 
commensurate value of Hi,, the vortex liquid (L) freezes 
first at Tg into the pinned smectic (S) state, followed by 
a second freezing transition at lower temperatures into 
the true vortex crystal (X). When He ^ 0, tilted smectic 
(TS) and crystal (TX) phases also exist. The TS-L and 
TX-TL transitions are XY-like, while the TS-S and TX- 
X phase boundaries are commensurate-incommensurate 
transitions (CITs) At larger tilts, the TX-TS and 

TS-L lines merge into a single first order melting line. As 
Hi, is changed, incommensurate smectic (IS) and crys- 
tal (IX) phases appear, again separated by CITs from 
the pinned phases, and an XY transition between the IS 
and L states. Note that the commensurate smectic or- 
der along the c axis is stable to point disorder because 
phonon excitations are massive (see below). This stabil- 
ity should increase the range of smectic behavior relative 
to the (unstable) crystalline phases when strong point 
disorder is present. 
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FIG. 1. Phase diagram for the pure vortex system as a 
function of temperature, T, and c-axis field. He- The fuzzy 
lines indicate first order transitions. A similar phase diagram 
holds in the 6Hb-T plane with the TX and TS phases replaced 
by IX and IS states, respectively. The merging of the IX-IS 
and IS-L lines need not occur in this case. 

We now proceed with the derivation of these results; 
further details will be given in Ref. The a, 6, and 

c, crystallographic axes are parameterized by the coordi- 
nates x,y, and z, respectively. The Cu02 layers are thus 
perpendicular to the z = c axis. The magnetic field is 
primarily along the y = b direction. When He = and 
Hi, is tuned to commensuration, de Gennes' freezing the- 
ory applies near Tj. Smectic order appears as a density 



wave, 

n(r)«noRe{l + $(r)e'9-}, (5) 

where uq is the background density, and q — 2TT/a is the 
wavevector of the smectic layering (with wavelength a). 
The complex translational order parameter <i>(r) is as- 
sumed to vary slowly in space. The superconductor is in- 
variant under translations and inversions in x and y, and 
has a discrete translational symmetry under z —^ z -\- s, 
where s is the Cu02 double-layer spacing. From Eqj^, 
these periodic translations correspond to the phase shifts 
$ — > $6"'". In the commensurate limit, a = ms, with m 
an arbitrary integer. The most general free energy con- 
sistent with these symmetries is 

^ = f l(V - ^A)$p + + 

_|($™ + $*™) + ...|, (6) 

where the coordinates have been rescaled to obtain an 
isotropic gradient term. The "vector potential" A repre- 
sents changes in the applied field 511 = SHi,y + HcZ, with 
Ax = 0, Ay = qHc/Hb, and = qSHb/Hb. The form 
of this coupling follows from the transformation proper- 
ties of $ . Additional interactions with long wave- 
length fluctuations in the density and tangent flelds 
are irrelevant to the critical behavior . 

When SH = A = 0, Eq.|| is the free energy of an 
XY model with an m-fold symmetry breaking term. A 
second order freezing transition occurs within Landau 
theory when v > and r oc T — Ts changes sign from 
positive (in the liquid) to negative (in the smectic). The 
renormalization group (RG) scaling dimension, A„i, of 
the symmetry breaking term is known experimentally in 
three dimensions as Am ~ 3 — 0.515m — 0Ab2m(rn — I) 
|]l5| . For J7i > rric ~ 3.41, the field g is irrelevant 
(Am < 0), and the transition is in the XY universality 
class. The magnetic fields used by Kwok et. al. 1^ cor- 
respond to m = 9 — 11 ||l^, well into this regime. The 
static critical behavior is characterized by the correla- 
tion length exponent v k, 0.671 ± 0.005 and anomalous 
dimension rj w 0.040 ± 0.003 |l|]. 

Deep in the ordered phase (r < 0), amplitude fiuctu- 
ations of $ are frozen out. Writing <i> = yJ\r\/ve^'"^^/°' , 
Eq.^ becomes, up to an additive constant, 

i^smectic = j d\\^^{Vu-Af -~gcos2T:u/sY (7) 

where k = 4,'K'^\r\K / a^v, g — g{\r\/v)"^/'^ , and the re- 
duced vector potential is .4 = A/q. The displacement 
field u describes the deviations of the smectic layers from 
their uniform state. The sine-Gordon term is an effec- 
tive periodic potential acting on these layers. As is well 
known from the study of the roughening transition JTsf , 
such a perturbation is always relevant in three dimen- 
sions. The smectic state is thus pinned at long distances 
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(i.e. the displacements u are localized in a single minima 
of the cosine) [ pj| . 

At lower temperatures, provided point disorder is neg- 
ligible, the vortices will order along the x axis as well 
^0| . This is once again a freezing transition at a sin- 
gle wavevector, and is described by a Landau theory like 
Eq.^ (with 5 = since the modulating effect of the un- 
derlying crystal lattice is much weaker in the x direction) , 
with XY critical behavior. 

Next, we consider the effects of (5H ^ 0. When the 
fields are small, A is an irrelevant operator in the smec- 
tic phase (as can be easily seen by replacing the peri- 
odic potential by a "mass" term cx v?). This implies 
that the smectic layers do not tilt under weak applied 
fields, i.e. d{dyu)/dH, 



-44 



d{B,)/dH, 



He 



Ha=o — 0. The full susceptibility 
is obtained from the expression 



B, = ^Im($*a„$) 
tlb 



-44,0 



(8) 



where 044,0 is the tilt modulus obtained from anisotropic 
GL theory (without accounting for the discreteness of 
the layers) and r" = d^r /dH^\He=o- Eq-H has a sim- 
ple physical interpretation. The first term is the con- 
tribution to Be from tilting of the layers (described by 
a phase rotation of $). Transverse field penetration at 
fixed layer orientation contributes via the second term. 
Such motion arises microscopically from a non-zero equi- 
librium concentration of vortices with large kinks ex- 
tending between neighboring smectic layers. Eq.H pre- 
dicts a non-divergent singularity C44(r) — C44(Ts) ~ 
\T — Ts\^~°' at the critical point, where a is the spe- 
cific heat exponent. At low temperatures in the smec- 
tic phase this crosses over to the much larger value 
C44 ~ [^/i±/UpkBT/{B(|)oms)]exp{E\\^/kBT), where the 
energy of a large kink is _Eik ~ yJe±Upms. 

As is well known from the study of the sine-Gordon 
model a larger incommensurability can be compen- 
sated for by energetically favorable "solitons" , or walls 
across which u — > u + s. Solitons begin to proliferate 
when their field energy per unit area CTficid ~ —k.As ex- 
ceeds their cost at zero field, ctq ^ (estimated from 
Eq.0). ^ ^ 

Physically, these solitons correspond to extra/missing 
flux line layers and walls of aligned "jogs" for SH along 
the b and c axes, respectively. In the former case, this 
leads to an incommensurate smectic (IS) phase, whose 
periodicity is no longer a simple multiple of s. For SU || z, 
the solitons induce an additional periodicity along the y 
axis. This tilted smectic (TS) phase has long range trans- 
lational order in two directions ||2^ . The analogous tilted 
crystal (TX) phase is qualitatively similar, but has long 
range order in 3 directions. 

Unlike the corresponding CIT in two-dimensional ad- 
sorbed monolayers | p^ , we find that energetic interac- 
tions between widely separated solitons dominate over 



entropic contributions ||ri|. The free energy density in 
the incommensurate phases is thus 



/e 



soliton 



(9) 



where a = Uficid + ctq < is the total areal free energy of 
the soliton; A and w set the energy and length scales of 
the soliton interactions. Minimizing Eq.^ gives a soliton 
separation I ^ wln{A/\a\) near the CIT. 

As the temperature is increased within the IS or TS 
phases, the system melts into the liquid. The IS-L and 
TS-L transitions are described by Eq.^ with g = (as 
can be seen by reversing the dilation discussed earlier and 
discarding oscillatory contributions to F proportional to 
g), and are thus XY-like. 

The shape of the CIT phase boundary is of particular 
experimental interest. In the mean field regime, this is 
obtained from the condition tr = as SH ^ |r|^, with 
Tmf = (™ ^ 2)/4. By the usual Ginzburg criterion, 



mean field theory breaks down for 



< 



{kBTv/K^/y. 



To determine the shape of the phase boundary in this 
critical regime, we follow the RG flows until \r\ is or- 
der one. Then SH^ - ^^^SH and gn ~ ^^"3, with 
^ ^ kl Rotational invariance at the rescaled fixed 
point {g — 0) implies that the field exponent is exactly 
Ah = 1 Using these renormalized quantities, we find 
Tcrit. = (I A™ I + 2)iy/2 « 4.9 - 7.2 for the fields used in 
Ref . 1^ . The IS-L and TX-L phase boundaries are non- 
singular and are determined locally by the smooth SH 
dependence of r. 

We next discuss transport measurements using a 
coarse-grained approach similar to flux line hydrodynam- 
ics |l4). The time evolution of the fluxon density n{r,t) 
is determined by the continuity equation 



dtn + • j„ = 0, 



(10) 



where V_l — {dx, dz), and is the vortex current, deter- 
mined by the constitutive equation, |14| 



5F „ 5F 



SF 



Ti^ = -nV^—+ndy—-TaV^—+nf, (11) 

on OT OTa 

where r is the coarse-grained tangent density, F is re- 
lated to the Bardeen-Stephen friction coefficient 7bs — 
noF, and the driving force is 



f = — J Ay 

c 



(12) 



Here J is the applied transport current density and //(r) 
is a random thermal noise. Upon projecting the critical 
smectic density modes near ±qz out of Eqjl^, we find 



7Bs9f$ = -4q' 



5Fr, 



(13) 
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where /j, = q(j)onQ/c and /^(k) = inoq'r]z{qz + k). Re- 
markably, Eq.|l^ has the same form as the model E dy- 
namics [ p3[ for the complex "superfluid" order parame- 
ter $, where now Jx plays the role of the "electric field" 
in the Josephson coupling. The actual electric field is 
~ jv,z'Po/c, leading via Eq.|l^ to 



no(t>Q 
' 2qc 



Im($*9t$) 



B 



+ (1 - 1*172) I 7^ ) Px.,„ 



(14) 



where Pxx,n is the normal state resistivity in the x direc- 
tion, whose appearance in the last term follows from the 
relation {nocjio / c)^ / Jbs ~ {B/Hc2)pxx,n- Eq|l| is inter- 
preted in close analogy with Eq.|^. The first term is the 
contribution to vortex flow from translation of the layers, 
while motion of equilibrium vortex kinks yields the sec- 
ond term. Such flow at "constant structure" is analogous 
to the permeation mode in smectic liquid crystals [|l2| . 

The presence of this defective motion implies a small 
but non-zero resistivity at the L-S transition. Near Ts, 
Eq.|l^ predicts a singular decrease of the form pxx (T) — 
Pxx{Ts) ^ \Ts - T|i"". At lower temperatures (but 
still within the S phase) transport occurs via two chan- 
nels. The permeation mode gives an exponentially small 
linear resistivity pxx ~ exp(— iJik/ZcB?") (above Tg, sin- 
gle layer kinks give pxx Gxpl—Ey^/k^T), with Ey_ « 
Ei]^/m). Non-linear transport occurs via thermally ac- 
tivated liberation of vortex droplets, inside which u (or 
Uz in the crystal phase) is shifted by s. Balancing the 
soliton energy on the boundary with the Lorentz en- 
ergy in the interior gives £ni ~ ^-^-^^l-^) ^ where Jc ~ 
(cj B^i^ngf'^'^islk^TY^'^ ■ A more detailed discussion of 
the full scaling form of £(J) will be given in Ref. 

In the TS phase, net vortex motion along the c axis 
occurs by sliding soliton walls along the h direction. 
The resulting electric field is proportional to J and the 
soliton density, leading to an additional contribution to 
the resistivity which vanishes at the CIT like p^*^''*"" ~ 
p^°^'*™/ln(A/|cr|). The situation in the IS phase is more 
complex, due to the possibility of a roughening transition 
for a single soliton wall, and will be deferred to Ref. . 

Lastly, we consider the effects of weak point disorder, 
which enters the free energy as a random field — 
J AnoVd(r)Re{$(r)e*''^}, where Vd{r) is a quenched ran- 
dom white noise potential. Such a perturbation alters the 
universality classes of the phase transitions in Fig.^, and 
renders all but the L and S phases glassy |]ll|. The sta- 
bility of the S phase to randomness (which takes the form 
Fd — J d'^^'^no^\r\/vVd{r) cos2n{u + z)/a in this regime) 
is a consequence of the phonon "mass" g/s'^. The nature 
of the glassy phases and altered critical behavior will be 
discussed in Ref. [ pl] |. 
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